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The electric microfield distribution at charged particles is studied for two-component 
electron-ion plasmas using molecular dynamics simulation and theoretical models. The parti- 
cles are treated within classical statistical mechanics using an electron-ion Coulomb potential 
regularized at distances less than the de Broglie length to take into account the quantum- 
diffraction effects. The potential-of-mean-force (PMF) approximation is deduced from a 
canonical ensemble formulation. The resulting probability density of the electric microfield 
satisfies exactly the second-moment sum rule without the use of adjustable parameters. The 
correlation functions between the charged radiator and the plasma ions and electrons are 
calculated using molecular dynamics simulations and the hypernetted-chain approximation 
for a two-component plasma, ft is shown that the agreement between the theoretical models 
for the microfield distributions and the simulations is quite good in general. 
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I. INTRODUCTION 

Because of the Stark effect, the fluctuating electric micro-fields created by the charged particles 
in a plasma influence its optical and thermodynamic properties. They affect the profiles of spectral 
lines (broadening and shift} and effectively lower the photoionization thresholds of atoms and ions 



immersed in a. plasma A comparison of experimental and theorel ica 



widths and shapes 



of the Stark-broadened spectral lines is widely used for plasma diagnostics \A al. 

P 

Under certain assumptions the observed spectral line shapes can be closely related to the 

electric microfield distribution at the radiating atom or ion (radiator) jf3,Q]. Within the quasistatic 
approximation the problem is then reduced to a determination of the probability distribution of 
the low-frequency component of the perturbing electric fields. This is mainly associated with the 
distribution of the heavier perturbing particles, i.e. the ions, whereas the electrons can be assumed 
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to adjust instantaneously to the configuration of the ions. 

Since the pioneering work of Holtsmark |g] , who completely neglected correlations between the 
particles (ideal plasma), many efforts have been concentrated on an improved statistical description 
of the microfield distribution. The first theory which goes beyond the Holtsmark limit and which 
is based on a cluster expansion similar to that of Ursell and Mayer |2| was developed by Baranger 
and Mozer j], ^| . In this approach the microfield distribution is represented as an expansion in 
terms of correlation functions which has been truncated on the level of the pair correlation. The 
latter is treated in the Debye-Hiickel form which corresponds to the first order of the expansion in 




the coupling parameter. The theory by Baranger and Mozer was improved by Hooper 
later by Tighe and Hooper Q|. Based on Broyles' collective-coordinate technique 
reformulated the expansion of the microfield distribution in terms of other functions by introducing 
a free parameter which was adjusted in such a way to arrive at a level where the resulting microfield 
distribution did not depend on the free parameter any more. A further improvement of this model 
was made in Ref. considering a Debye-chain cluster expansion. Afterwards the Baranger- 
Mozer second order theory was extended by including higher order corrections, like the triple 
correlation contribution jl?! . Iisl |. However, it was argued that such a method is only valid for 
low-density, high-temperature plasmas, i.e. at small coupling parameters, where the correction to 
the Holtsmark distribution, corresponding to the first term in the series, is small. In the limit 
of very strong coupling and without screening Meyer's harmonic oscillator model is applicable 
jigl |. in which every ion is assumed to oscillate independently of the others around its equilibrium 
position at the ion-sphere center. The first theory capable to provide reliable numerical results 



for strongly coupled plasmas, known as adjustable- 



was proposed by Iglesias, Lebowitz et al. 
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parameter exponential approximation (APEX), 
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This phenomenological but highly 



successful approximation is based on a special parameterization of the electric microfield produced 
on a radiator. It involves a non-interacting quasiparticle representation of the electron-screened 
ions, designed to yield the correct second moment of the microfield distribution. APEX was first 
developed for three- and two-dimensional Coulomb systems 12 ll 



Coulomb systems and ion mixtures 2Jj 



2J. (See also Ref. 



, |22j and later adapted to screened 
for the corrected version of APEX 



for a neutral radiator). Another approach providing reliable numerical results for the strongly 
coupled plasmas was proposed by Iglesias 2f|. Following the idea of Morita 27] on the similarity 
of the representation of the microfield distribution to that of the excess chemical potential, Iglesias 
reduced the problem to a determination of the radial distribution function (RDF) for a fictitious 
system with an imaginary part in the interaction energy. Employing this idea Lado and Dufty 



28 



29, 30] developed an integral equation technique for calculating the RDF and good agreement 



was found with computer simulations. It is now possible to calculate the microfield distribution 

la5mfl 0QQQQ. 

These methods allows to study the effects of microfield nonuniformity 



from Monte Carlo (MC) or molecular-dynamics (MD) simulations of p 



properties of the electric microfiek 



distribution in electron plasmas 25 



36, 371 and the dynamical 

n " 

38] as well as to simulate the high-frequency microfield 
]. With these powerful tools one can check the accuracy of 



theoretical 
(see, e.g 



models and establish asymptotic or analytic fitting formulas suitable for applications 



40, 



4l| and references therein). 
Until now most work was done on either electronic or ionic one-component plasmas (OCP) 
neglecting the influence of the attractive interactions between electrons and ions. Here we treat ions 
and electrons on an equal footing by concentrating on two-component plasmas (TCP). Previously 
this has been done in Ref. 3] for partially degenerate electrons. In particular, the low-frequency 
component of the microfield was calculated within the linear response treatment taking strong 
correlations into account via local field corrections. Also the problem of attractive interaction has 
been considered for single but highly charged impurity ion immersed in an electronic OCP (see, 

n 

e.g. Ref. |43j for a recent review of these cases). 

In the present paper we study strongly coupled systems, i.e. a highly charged radiator in a TCP 



of classical 
As in Ref. 



nondegenerated) and strongly correlated particles beyond a perturbative treatment. 



421 ] the presented theoretical scheme is based on the potential-of-mean-force (PMF) 



approximation which exactly satisfies the sum-rule requirement arising from the second moment of 
the microfield distribution without introducing adjustable parameters. Another important ingre- 
dient is the electron-ion attractive interaction which drastically changes the physical properties of 
the system as compared to classical OCPs (see, e.g., ^jj]). This may cause significant changes in 
the microfield distribution on either neutral or charged radiators. But the thermodynamic stability 
of a TCP requires some quantum features for the electron-ion interaction at short distances. Here 
we focus on an application of classical statistical mechanics and MD simulations which is enabled 



by using a regularized ion-electron potential where the divergence at the origin is removed 



see also [lil 1 1 ■ For a review. 

The paper is organized as follows. In Sec. II, we define the systems and parameters of inter- 
est as well as the theoretical model to calculate the microfield distribution in a TCP. The exact 
second moment for the charged radiator is calculated in Sec. III. The theoretical schemes applied 
previously to either electronic or ionic OCPs are generalized to TCPs in Sec. IV. In particular, 
we consider the Holtsmark distribution, express the microfield distribution through the pair dis- 



27, 



44 



tribution functions, and deduce the PMF approximation from the classical canonical ensemble. 
Furthermore we construct a theoretical approach based on the exponential approximation where 
the effective electric fields are calculated on the basis of the PMF approximation and the pair corre- 
lation functions. In Sec. V we consider the hypernetted-chain (HNC) integral equations technique 
to calculate these functions in a two component plasma. In order to test the theoretical models 
we carried out classical MD simulations to calculate both the pair correlation functions and the 
microfield distribution. Technical aspects and the numerical results are presented in Sec. V. These 
results are summarized in Sec. VI. Some details of the calculations are described in the Appendix. 



II. MICROFIELD DISTRIBUTION IN A TCP: THEORETICAL BACKGROUND 

A. Basic parameters for the TCP 

We consider a neutral and isotropic two component electron-ion plasma consisting of iVj ions 
and iV e electrons at a temperature T in a volume The particles are assumed to be classical 
and pointlike. The average densities, charges and masses of the ions and electrons are rij = Ni/O,, 
n e = N e /Q, and Ze, —e and rrii, m, respectively. We assume that the density of radiator ions 
are small, tir <C nj ;e and thus consider only one radiator ion with charge Z^e in our calculations 
(throughout this paper the index R refers to the radiators). Because of the charge neutrality we 
have NiZ — N e + Zr = 0. In the thermodynamic limit (iVj ;e — > oo and Q — > oo) this is equivalent 
to n e = riiZ. 

We now introduce the Coulomb coupling parameters T a p which play an important role for 
characterizing the properties of a TCP. Introducing the Wigner-Seitz radii, i.e. the mean electron- 
electron, electron-ion and ion-ion distances through the relations, a~ 3 = 47rn e /3, a~ 3 = 47rn/3 and 
a~ 3 = 47rnj/3 (where n = n e + rii is the plasma total density) these parameters are defined as 

e| Ze 2 Z 2 e 2 

-Tee — J TFii ^ei — J T^i ^ii — ; ™j (1) 

a e KB-l aksl CLikBi 

respectively, where e| = e 2 /^tteq. Note that 

ee " [z 2 {z + i)f u ~(z + if r { ) 

In a hydrogen plasma with Z = 1 W6 obtain r ee — — 2- 1 / 3 T ei while in a plasma with highly 
charged ions [Z » 1) T lt = Z 2 l 3 T ei and r ee = T e i/Z. For Z ^ 2 the coupling parameters satisfy 
the inequality T ee < T e i < Tu. 



Here we consider the pair interaction potential e^q a qpu a p (r) with a; ft = e,i,R, q e = — 1, 
qi = Z, q R = Z R , and 

u aP {r) = - r (l-e- r l 5 ^) (3) 

which is regularized at small distances due to quantum-diffraction effects. In this paper we assume 
that the Coulomb potential is cutoff at the thermal de Broglie wavelengths, 5 a p = (h 2 / ' /J-apkBT) 1 ^ 2 
where [i a p is the reduced mass of the particles a and (3. For large distances r > 5 a p the potential 
becomes Coulomb, while for r < 5 a p the Coulomb singularity is removed and u^ih(Q) = ll8 n H. By 
this the short range effects based on the uncertainty principle are included 27|, 143 . |4J, |45 1 . 

For a classical description of a plasma the electron degeneracy parameter O e , i.e. the ratio of 
the thermal energy and the Fermi energy must fulfill e = ksT/Ep > 1. Or, alternatively, the 
electron thermal wavelength should be smaller than the electron-electron mean distance, S ee < 
2(4/9vr) 1/3 a e ~ 1.04a e . Since an ion is much heavier than an electron this condition is usually 
fulfilled for ions. We note that 5u <C <5 e i and 5 ee ~ 2 1 / 2 5 e i since \x e i ~ m. Therefore one can expect 
that the regularization given by Eq. (3) is less important for ions than for electrons. Furthermore, 
scattering of any two particles is classical for impact parameters that are large compared to the 
de Broglie wavelengths. Typical impact parameters are given by the Landau lengths, Xia/3 = 
e| If/o^l /ksT. Its ratio to the de Broglie wavelengths is given by 

^Laf3 ^ \q a qp\ a e 2 s \q a qB\u a 8(0) 
Oot/3 Z d a[3 k B T 

This is also the maximum value of the interparticle interaction energy in the units of fcgT, where 
a ee < a e i <^ an and a e i ~ 2 1 l 2 Za ee . Classical description of the scattering events in the TCP is 
valid if c ee > 1. This can be alternatively written in the explicit form k R T < lRy Combining 
this condition with the one considered above we finally obtain the temperature domain where 
the classical treatment is adequate, Ep < k R T < 1 Ry. This condition occurs at lower densities 
of electrons. Since the parameter 5 ee increases with electron- ion Coulomb coupling the classical 
condition a ee > 1 implies that the state with stronger T e i behaves more classical as discussed in 
Ref. Q- 

B. Microfield distribution formulation within thermodynamic canonical ensemble 



The electric microfield distribution (MFD) Q (e) is defined as the probability density of finding 
a field E = e at a charge Z R e, located at ro, in a TCP with iV, ions and N e electrons. This system 



is described by classical statistical mechanics in a canonical ensemble of (iVj + N e + 1) particles, 
and temperature T. The normalized probability density of the microfield s in the thermodynamic 
limit is then given by 

w Jn 

where 0t = 1/ksT, and T e = {ri, r2...rjv e }, % = {Ri, R2...RjVj} are the coordinates of electrons 
and ions, respectively. Here 

W = I e-^ u ^' T ^dr dT e dTi (6) 
Jn 

is the canonical partition function and U (%,,%, tq) is the potential energy of the configuration 

U (%,%, r ) = U ee {%) + U u (%) + U ei (T c , %) + U eR {%, r ) + U iR (%, r ) (7) 

with electron-electron, ion-ion, electron-ion, electron-radiator and ion-radiator interaction terms, 
respectively. Assuming spherical symmetric interactions between the particles the interaction terms 
in Eq. (7) can be represented as 

U a p {%,T p ) = # af} q a qpe 2 s £ u a(i ^ - rf } |) , (8) 

a,b 

U aR (T a ,r )=q a Z R e 2 s Y, u *R(\ro-ri a) \) (9) 

a 

in terms of the pair interaction potentials u a p (r) and u aR (r), where a; (3 = e; i, $ ee = •da = 1/2, 
$ e i = 1, = r a , r$ = R a . In Eq. (8) the sum is restricted to a ^ b for like particles, a = f3. The 
total electrical field E (T e , %,ro) acting on the radiator is given by the superposition of electronic 
and ionic single-particle fields 

B(%,T t ,r ) = -—V U = E e (%,r ) + E i (T i ,r ) (10) 
Z R e 

with 

N a 

E a (T a ,r ) = ^E Q (r -ri a )). (11) 

a=l 

As E e (r) = j,E e (r), Ej (r) = ^Ei (r), we obtain for the electronic and ionic single-particle fields 
E e (r) = eFu' eR (r), E{ (r) = —Zepu' iR (r), where the prime indicates derivative with respect to r, 
and eF = e/AneQ. 



The spherical symmetric interaction between plasma particles allows to introduce the normalized 
microfield distribution P(e) = 4ire 2 Q(e). It is useful to consider the Fourier transform of Q (e) 
defined by 

T(K) = j Q (e) e iK - £ de = <e iK - E > . (12) 

Here (...) denotes a statistical average. Again we note that due to the isotropy of the system the 
Fourier transform of the MFD must behave as 

.2 /"oo 



T(K)= P (e) j (Ke) de, P (e) = — / T (K) j (Ke) K 2 dK, 

J0 7T Jo 



(13) 



where jo (x) = sin x/x is the spherical Bessel function of order zero. The coefficients of the expansion 
of the function T(K) at K — > yield the even moments of the microfield distribution, 

^) = l-^<£ 2 > + ^<£ 4 )-- (14) 
The similar expansion for the function C (K) defined by T(K) = e~ c ^ yields 

+ ... (15) 



k 2 ^ 4 r,„ 2 \2 3 



Therefore the Fourier transform of the MFD can be interpreted as a generating function for mi- 
crofield even moments. Moreover, Eqs. (14) and (15) suggest a simple criterion for the existence of 
even moments. In particular, the second moment of the MFD exists if the function C (K) and its 
first and second derivatives are regular at the origin. Eqs. (5)-(15) then describe the total MFD at 
the position ro of the radiator generated by both the statistically distributed ions and electrons of 
the TCP. Since we are interested to calculate the MFD, Eq. (5), in an infinite system the statistical 
average of any quantity becomes translationally invariant with respect to ro, and the location of 
the test charge may be taken as the origin without loss of generality. 



III. SECOND MOMENT 



A knowledge of moment sum rules is often useful in developing approximation schemes for 
fluids and plasmas. The moments of the MFD fix the shape of the distribution and involves 
some useful information about the system. For example, the exact second moment has been 
previously incorporated into the calculation of the MFDs in the APEX scheme. Here, we derive 
exact expressions for the second moment of the MFD on charged radiators. Note that in general 
the existence of the second moment requires that the MFD decays at large electric fields faster 
than e -3 . 



Let us consider the exact expression for the second moment of the microfield distribution in the 
TCP and for a charged radiator. The second moment may be written in the form 

<£2> = ^( (v ° c/)2 )' (16) 

where Vo is the gradient with respect to ro and the average is over the canonical ensemble defined 
in Eq. (5). Noting that e"^ TU (V U) = -k B T (V e~ l3 ' rU ), substituting this relation into Eq. (16), 
integrating by parts, and setting the surface terms equal to zero yields 

<*">-(i^< v * 7 >--^« v '- B »- (17) 

We now use Eqs. (10) and (11), the relation V • E„(r) = ((foe^/r 2 ) u a (r), where u a (r) = 
— [ r2u ' a R ( r )]'i an d translational symmetry. This yields 

( E 2\ = k BTn e 
V 1 Z R e 

The functions g a n{r) are the pair correlation functions between radiator and the plasma particles, 
where n a g a R(r) is the density of plasma particles a at a distance r from the radiator. These 
functions can be represented as 

9eR{n) = %! e-WKMdlWdTi, (19) 



oo poo 

u e (r) g e R (r) dr - ui (r) g iR (r) dr 
uo Jo 



(18) 



9iR (Ri 



f e -PTU(T e ,7l) dTdT (l) 



(20) 



Here dT^ = Ilo=s+i ^ r «" ; ^ s ^ ne reduced volume element in a phase space which does not involve 
the particles 1,2, ...s of plasma species a. The interaction potential energy, U (%,%), does not 
depend on ro- The pair correlation functions given by Eqs. (19) and (20) describe the coupling 
between radiator ion and plasma particles. For a vanishing radiator-plasma coupling, e.g. for a 
neutral radiator the pair correlation functions behave like g a R — > 1. If the radiator is a particle of 
plasma species (3 these correlations functions coincide with the radial distribution functions (RDF) 
of bulk plasma, g aR = g a/3 . 

The second moment for the regularized Coulomb interaction (see Eq. (3)) is with u a (r) = 



-r/5 aR 



k B Tn e 
ZrSq 



- r/SeR geR (r)rdr 



%RJ 



- r/SiR 9iR {r)rdr 



(21) 



Usin g a bare Coulomb interaction u a (r) = 5 (r) in Eq. (18) one recovers the result obtained in 



Ref. 



42] 



^ = iSt [geR (o) " 9iR (o)] ' (22) 

which can also be obtained from Eq. (21) by taking the limits 5 e R — > 0, Jjjj — > 0. For Zr > 0, we 
may assume that gm (0) = if quantum-diffraction effects are negligible for the ions, while g e R (r) 
diverges at small distances for a bare Coulomb potential. This indicates that the second moment 
of the microfield distribution does not exist for a classical Coulomb TCP. But in the OCP limit 
QeR (0) = 1 one recovers the result (E 2 ) OCp = kBTZrii/ZjiEo for the classical (ionic) OCP j^. 



IV. APPROXIMATE CALCULATIONS OF THE MFD 



In this section we generalize the existing theoretical approaches developed originally for a OCP 
to a two component electron-ion plasma. For practical applications we will consider the exponential 
approximation considered in Ref. ,46;], and, as a simple but useful example the Holtsmark limit for 
the MFD in a TCP. 



A. Ideal plasmas: Holtsmark distribution 

We first consider the microfield distribution in an ideal TCP with T ee ,T e i,Ta — > 0, i.e. in the 
high temperature regime T — > oo. In this case Eq. (12) yields 

T {K) = J! jl - ^ jf [1 - jo (KE a (r))] r 2 drp . (23) 

In the thermodynamic limit (N a , O, — > oo, N a /Q, = n a = const), and recalling that T{K) = e~ £ ^- > 
we obtain from Eq. (23) 

/•oo 

C(K)=Y / 4vm Q / [1 - j (KE a (r))] r 2 dr. (24) 

We study this expression for two types of interaction potentials. 

(i) For a bare Coulomb interaction Eq. (24) yields C (K) = (KEh) 3 ^ 2 , where Eh is the Holts- 
mark field for a TCP, E% 2 = E^+E^f. Here E He and E m are the electronic and ionic Holtsmark 
fields, respectively, Ejj e = Cep/a 2 , Eni = CZep/af with C = (87r/25) 1//3 . Since Eu e = Z~^ 3 Ejji 
the electronic and ionic components of a hydrogen TCP contribute equally to the Holtsmark field. 
For a completely ionized TCP with highly charged ions the ions dominate Eh- The definition of 



the Holtsmark field Eh for a TCP is equivalent to the obvious relation n = n e + rij and can be 
represented as 

2/3 

(25) 

with an effective charge Z. For a hydrogen TCP with Z = 1 also Z = 1. In other cases the effective 
charge increases with Z and behaves as Z = Z 1 / 3 for large Z. Thus the ideal two-component plasma 
can be regarded as an ionic OCP with effective ionic charge Z. 

Since the function C(K) has a singularity at K = it cannot be expanded there and the second 
moment does not exist. The microfield distribution is given by Pjj (E) = H (rj) /Eh in terms of 
Holtsmark's function H(r]) 

2ri f°° _ 3/2 

H (v) = — / e ~ x srn ( 7 ? x ) x d x (26) 

with Tj = E/Eh- Note that the Holtsmark distribution for the TCP has the same functional form 
as either the ionic or the electronic OCP. The only difference is the definition of the Holtsmark 
field. Since the electronic or ionic Holtsmark fields may significantly differ from Eh the shape of 
the MFD for a OCP and a TCP may strongly differ from each other even for ideal plasmas. 

(ii) For the regularized Coulomb interaction given by Eq. (3), C{K) (from Eq. (24)) and all its 
derivatives are regular at K = 0. Hence, all moments of the microfield distribution exist. This 
indicates that for large electric fields the microfield distribution must decay exponentially. The 
second moment can be obtained from Eq. (24) if we recall that for K -> 0, C(K) ~ (if 2 / 6 ) (£ 2 ); 
thus 

<£ 2 > = 2,n e el i^ + iPj. (27) 

\0 e R OiRj 

For large electric fields the main contribution to the microfield distribution comes from small K 
and we obtain the asymptotic behavior 

where (-E" 2 ) is given by Eq. (27). For large K (K — ► oo) the function C(K) (Eq. (24)) behaves as 
for the bare Coulomb interaction C{K) ~ (KEh) 3 ^ 2 ■ Hence, the microfield distributions for the 
ideal plasmas with bare and regularized Coulomb potentials behave similar at small electric fields. 



1/3 



Z(l + Z 1 / 2 ) 
Z + l 



B. Expression of the MFD through pair functions 

It was first noted by Morita j^] that the virial expansion of the Fourier transform of the MFD 
T(K) is formally similar to that of the excess chemical potential. This was previously used to 



express T(K) in terms of an effective RDFs (see, e.g., |26j) involving the radiator and one of the 
plasma particles. To generalize this method to the TCP we follow the procedure [6-30] and consider 
the logarithmic derivative of Eq. (12) 



aCGO ( ( K ■ E J e iK ' E 



dK 



iK-E\ 



iK-Y,n a J drB a (r) [G aR (r, K) - 1] . (29) 



Here E e (r) and Ej (r) are the single-particle electronic and ionic electrical fields introduced above, 
and K is a unit vector in the direction of K. G e R ( r > K) and GiR (r, K) represent the pair correlation 
functions between the radiator and the plasma particles in a fictitious system whose interaction 
potential is given by the complex quantity U (T e , %, K) = U(T e , %) — i (ksT) (K • E), i.e., 

G eR (r 1: K) = ^ J^e-^^drPdT,, (30) 

GiR (Ri, K) = y^j n e-^%,K )d%dr m (31) 

with the generalized, reduced partition function W(K) = W(K) = W(e lK ' E ). It can be easily 
checked that this function is real. In general these correlation functions are complex and satisfy 
the symmetry relations G a R{— r, K) = G^ R (r,~K) and Q a R (— r, — K) = C? Q jj(r,K), where the 
asterix denotes the complex conjugate. The correlation functions in the fictitious system are 
not spherical symmetric. At K = they coincide with g a R(r) given by Eqs. (19) and (20). 
The complex correlation functions Q a R can be expressed through two functions (r, K) and 
4 0) (r,K) = i£i 0) (r,K), 

Q aR (r, K) = (r, K) exp \iK ■ £f (r, K)] , (32) 

where G^r ( r ' K) an< ^ ( r ' are spherical symmetric real functions. With this choice the 
functions G a R ( r > K) automatically satisfy the symmetry relations. In the limit K — > 0, we have 
also ^°](r,0) = 9aR(r). Inserting Eq. (32) into Eq. (29), integrating over K and taking into 
account that £(0) = we obtain 

poo f-K 

C (K) = 4vr £ n Q ^ ^ a (r) r 2 dr ^ ^ (r, A) ^ (a4°) (r, A) J dX, (33) 

where ji (x) = —j' (x). Eq. (33) is an exact result which allows to express the MFD through 
complex pair correlation functions (or, alternatively through two real functions). In addition 
Eq. (33) yields the exact second moment given by Eq. (18). 

The problem is now the evaluation of these correlation functions. Eq. (33) requires that the 
complex correlation functions has to be known in the interval from to K. One possibility is to 



apply the integral equation technique with the complex interaction energy introduced above. Such 
an approach has been previously employed for a OCP [2^ . I2H I30I and shows good agreement with 
computer simulations. Here we adopt the exponential approximation (see, e.g., 



25, 



26, 



46 



theory 



21 



, 22, 



;11L Willi 



and generalize it to the TCP. This method is based on thermodynamic perturbation 
471 ] . The system with the potential energy IA{K = 0) = U is chosen as reference system and 
its structure is assumed to be known to a good approximation. The perturbation potential is then 
given by U\ = —i (k R T) (K • E) and we expand the correlation functions, Eqs. (30) and (31), with 
respect to U\. Within first order we obtain Q aR (r, K) ~ g aR (r) [1 + iK • £ a (r)]. Here g aR (r) are 
the actual RDF in the real system, Eqs. (19) and (20), and £ a (r) = ( r >0). The electric fields 
£ a (r) may be interpreted as effective electric fields in the fictitious system which are independent 
of K. Taking into account that (E) = we obtain 

S a (r) = E a (r) + - J— Y^np [ d^Ep (n) [g a/3 (|r - n|) - 1] . (34) 
9aR (r) ^ J 

Comparing Eq. (34) with Eq. (3.8) of Ref. [sjj for the case of a OCP we remark that our present 
derivation yields an additional factor \/g aR (r) in front of the second term. Since g a fs( r ) depend 
only on |ri — r2 1 the effective electric fields in Eq. (34) can be represented as £ a (r) = r£ a (r). The 
S a (r) can be expressed by the pair correlation functions and the single-particle potentials u a R (r) 
(see Appendix A for details). Alternatively the Fourier transformed single-particle electric fields 
can be written as E Q (k) = \iE a (k) which allows to express the effective fields through the static 
structure factors S a p (k). 
We make now the ansatz, 

G aR (r, K) = g aR (r) exp [*K • £ a (r)] (35) 

and then integrate Eq. (33) with respect to A, to find 

C (K ) = V 47m Q [°° E a (r) 1 ~ j °} K . £ . a {r)) g aR (r) r 2 dr. (36) 
a Jo t-a \r) 

The second moment within the exponential approximation can be found from Eq. (36) at K — > 
and results in 

/>oo 

(E 2 ) = V 4vrn a / E a (r) £ a (r) g aR (r) r 2 dr. (37) 

This must fulfill the exact second moment of the MFD given by Eq. (18) or (33), which is not 
affected by either the assumption (35) or its first order Taylor expansion with respect to K. 



The APEX approach was originally developed for the classical ionic OCP with bare Coulomb 
interaction. In order to fulfill the exact second moment (-E 2 ) OCp = kBTZrii/Zneo, Eq. (37) must 
take the form 



r 8{r)g R {r)dr = ^. (38) 



In Ref. 



2l| the effective field 8 (r) is assumed as a Debye-Hiickel like screened interaction with 
unknown screening length. This free parameter is then adjusted in such a way to satisfy Eq. (38). 
The resulting predictions of APEX for the probability densities show excellent agreement with 
numerical simulation data for the OCP. However, difficulties appear when one attempts to extend 
the APEX scheme to a TCP, e.g. by assuming a Debye-Hiickel like interaction separately for 
the electrons and the ions and introducing two adjustable screening lengths. Then the sum rule 
Eq. (38) with the exact second moment {E 2 ^ becomes ambiguous as it allows for many different 
choices of the adjustable screening lengths. This can be cured for ionic mixtures by demanding 



that the second moment rule is satisfied species by species (see, e.g., |23j, |24|). But this cannot be 
employed for a TCP with attractive electron-ion interactions. Here the Debye-Hiickel ansatz for 
the electronic effective field is physically incompatible with Eq. (37) as discussed in Ref. j^ . 
We instead apply the potential of mean force (PMF) approximation 42 , |3] which expresses 



the effective electric fields through the logarithmic derivative of pair correlation functions 

£ a (r) = ^-^- [In g aR (r)]. (39) 

Introducing Eqs. (39) in Eq. (37) automatically satisfies the sum-rule (18) without any adjustable 
parameter. Relations (39) can be deduced from Eqs. (19), (20) and (34). To show this we con- 
sider the pair correlation functions given by Eqs. (19) and (20). It is clear that in the thermo- 
dynamic limit these expressions are translationally invariant with respect to — > + 
(a = 2, 3, N a ). Making these transformations and calculating the logarithmic derivatives of 
the pair correlation functions yield Eq. (39), where the effective fields are given by Eq. (34). In 
addition, Eq. (39) can be interpreted as a integro-differential equation for determining the pair 
correlation functions, g a R- Thus, if the g a R( r ) are known the MFD with the exact second moment 
can be calculated using Eqs. (13), (36) and (39). This approach based on the exponential and the 
PMF approximations is abbreviated as PMFEX in the following. 

We summarize this section by the following remarks. The possibilities of the PMF approxima- 
tion have already been noted by Alastuey et al. [22J. They found a superiority of the APEX to 
the PMF approximation since the former reproduces the simulation data for classical ionic OCP 



more accurately than the latter. We have confirmed this by own investigations on the OCP. For 
the TCP the outlined PMFEX approximation agrees quite well with the MD simulation results, 
as we will show in the next section. 



V. RESULTS 



In Sec. IV we introduced and outlined the PMFEX approximation which links the MFD to the 
RDFs. To obtain explicit results for the MFD the corresponding RDFs have to be determined first. 
This will be done by solving numerically the Hyper-Netted-Chain (HNC) integral equations for the 
TCPs under consideration. The HNC method and the PMFEX approximation are tested both by 
comparison of the resulting RDFs and MFD with those obtained by classical MD simulations. We 
have done that for a wide range of coupling parameters T e i and for two specific rather distinct cases 
H + (n e = rtj) and Al 13+ (n e = 13nj) TCPs with symmetric and asymmetric density distributions 
between plasma species, respectively. For simplicity we assume bare Coulomb electron-electron 
and ion- ion interactions with 5 ee ~ and 5a ~ while the parameter 5 e i/a = 5 scaled in the 
Wigner-Seitz radius a = [Air (n e + n,) /3] _1//3 varies from 0.1 to 0.4. 



A. Numerical treatments 

To determine the RDFs g a R{f) the HNC equations (see, e.g., Refs. |4Sl l49| ) 

1 + Kp (r) = exp [h a p (r) - c a p (r) - PTq a qpe 2 s u a p(r)] , (40) 
and the Ornstein-Zernike equations 

K/3 (r) = c a p (r) + y^ j n a dr'c aa (|r - r'|) h a p (r) (41) 

a J 

for the total correlation functions h a p (r) = g a p (r) — 1 and the direct correlation functions c a p{r) 
are considered. This has to be done for a three-component system of electrons, ions and the radiator 
in general. Here we assume that the radiator is one of the plasma ions Zr = Z, i.e. g a n(r) — > g a i(r), 
which reduces Eqs. (40) and (41) to the HNC-scheme for a TCP with mutual interactions u a p(r) 



ly_by an iterative 
Within our 



(see Eq. (3)). The resulting coupled equations (40) and (41) are solved numerica 
scheme which closely follows the implementation discussed in detail in Ref, 
numerical treatment a parameter regime with a e i = a = T e i/8 < a c (Z,5) is accessible, where the 
critical value a c for 5 = 0.1, 0.2 and 0.4 takes the values a c ~ 8.32, 8.5, 13.4 and a c ~ 7.33, 6.66, 7.0 



for H + and Al 13+ TCPs, respectively. Beyond this value the numerical procedure does either not 
converge or ends up in unphysical solutions. A similar behavior has been reported in ^| for the 
case of an ion embedded in electrons. With the RDFs provided by the HNC-scheme the MFD, 
i.e. P(E), is then calculated via Eqs. (13), (36) and (39) by standard numerical differentiation and 
integration methods |5ll |. 

In the MD simulations the classical equations of motion are numerically integrated for Ni ions 
and N e electrons interacting via u a p(r) and contained in a cubic cell with periodic boundary 
conditions. To account for the long range of the Coulomb interaction the forces are calculated 



by an Ewald sum 
Verlet algorithm 



52 



54, 



53j. The numerical propagation is accomplished by a standard Velocity- 



5 a] extended by a hierarchical treatment of close colliding particles which 



are propagated as subsystems (see [56j for details) and using an adaptive time step. Such MD 



simulations have already been extensively tested and successfu 



57, 



l y a pplied for investigations of the 



581. 



dynamic response of a TCP with regularized potentials, see 

The actual simulations run with N = Ni + N e = 2002 particles and proceed in two phases. An 
initial equilibration starts from a random sampling of positions and velocities and relaxes towards 
the equilibrium distribution of desired temperature by dynamic propagation with velocity rescaling. 
The subsequent simulations are performed in the microcanonical ensemble, where their accuracy 
and stability can be monitored using the total energy. The MFD and the RDFs are sampled during 
the simulations from the known forces on the particles and their positions as a time average over 
the total running time r which was typically r ~ 700u~ le , where oj p \ tC = (n e e 2 /meo) 1//2 is the 
electronic plasma frequency. 

By the MD simulations basically all correlations and many-body effects of classical many-body 
systems can be taken into account. Limitations arise mainly from the finite particle number and 
the system size, e.g. in connection with the screening of the interactions on a typical screening 

1/2 

length Xd. Since Xd should be smaller than the size L of the simulation box and L/Xd oc T ei , 
the MD-technique works here more favorable at large coupling (r e j > 1) while the limit of weak 
coupling (r e j <C 1) requires a strong increase of the simulation box, i.e. of the particle number. 



B. Correlation functions 



In Figs. 1-3 we compare the RDFs calculated either from the HNC scheme or MD simulations 
for a H + -plasma. Only g e i(r) and gu(r) are plotted as g e e{f) = 9n(r) for hydrogen both in the 
HNC and MD treatment (within numerical fluctuations). In Fig. 1 we explore the dependence in 



the regularization parameter 5 at a fixed electron-electron coupling T ee = 0.1, while in Fig. 2 T ee is 
varied. Both approaches agree perfectly in the range of parameters covered in these figures. Due 
to the regularization of the ion-electron interaction the RDF g e i(r) is finite in the limit r — > 0. A 
non- linear Debye-Hiickel approximation for g e i(r) has been proposed in Ref. |43j. Adopting this 
estimate for TCPs we obtain g e i(0) — exp (T e i/lZ5) with 1Z = 1 + 5 {ZT e i) 1 ^ 2 , where the dependence 
on the ion charge Z is included in the coupling parameter r e j. The RDFs show indeed the expected 
growth of correlations with increased coupling and decreased regularization parameter. For very 
strong electron-electron coupling deviations between the HNC scheme and the MD simulations 
begin to appear as shown in Fig. 3 for T ee = 4.0 and o"h+ — 12.5. There we are at the edge of the 
HNC convergence region and classically bound states show up in the MD simulations. 

The symmetry between the correlation functions g e e( r ) an d gu(r) breaks down for Al 13+ -plasma. 
In Fig. 4 we compare the HNC and MD radial distribution functions for fixed T ee = 0.01 varying 
the regularization parameter 5. Since Tu ^> T ee strong correlation effects are expected for ga. The 
increasing "correlation hole" is clearly visible in Fig. 5 where the HNC and MD radial distribution 
functions ga(r) and <7 e i(r) are plotted for a fixed 5 = 0.4 and varying coupling strengths T ee . Again, 
for all these parameters the HNC scheme agrees perfectly with the MD simulations. As shown in 
Fig. 6 deviations occur in the electron-electron RDF g e e( r ) at small r for strong coupling T ee = 0.2, 
<r A1 i3+ ~ 6.6. These are due to the enhancement of the electronic density around an ion, which also 
increases the probability of close electronic distances and results in the maxima in g ee at distances 
r < a. This effect is obviously overestimated in the HNC approach and it is more pronounced for 
highly charged ions like Al 13+ and less important for H + . The regularization of the electron- ion 
interaction has no visible influence on the correlation functions g e e( r ) and gu(r) (see Figs. 1 and 4). 

C. Microfleld distribution 

We now turn to the MFDs at the charged reference point which is chosen to be one of the 
plasma ions, Zr = Z. For our analysis it is instructive to consider first the second moment which 
can be used to check and compare the different treatments PMFEX, HNC, MD, and can provide 
some information about the shape of the MFD, although this is not sufficient to construct it. For 
a bare ion-ion Coulomb interaction the second term in the rhs of Eq. (21) vanishes and the second 
moment (E 2 ~) receives contribution only from the first term involving g e i. In the limit of an ideal 
TCP (r a/ 3 — > 0) g e i(r) can be replaced by unity. This yields (E 2 ) = ksTrie/ Zeq = (3 /Ta) = 
(3/Z 5//3 r ee ) with Eoi = Zep/af which is similar to the second moment obtained for the ionic 



OCP (see, e.g., (21]]). In this sense the ideal TCP behaves like an ionic OCP with F = Fa. The 
second moments calculated from Eq. (21) using a HNC radial distribution function g e i(r) are shown 
in Figs. 7 and 8 as a function of F ee for hydrogen and aluminum TCPs, respectively. The dashed 
straight lines represent (-E ,2 ) for the ideal system. The other curves are calculated for different 
5 up to the critical values a c = F^ /S = (f^J /6j [Z 2 (Z + l)) 1 ^ 3 introduced above. For small 



r ee the deviations of the second moment from {E 2 ^ Q are small and the second moment decreases 
approximately as 1/T ee . But unlike (-£/ 2 ) of an ionic OCP, it increases again with the coupling 
parameter F ee due to the strong attractive ion-electron interactions. 

The normalized MFDs from PMFEX and MD are compared in Figs. 9-18 where the electric 
microfields are scaled in units of the Holtsmark field Eh (see Eq. (25)). For each distribution we 
have also calculated the second moment as a control parameter and found a quite good agreement 
between Eq. (21) and the MD simulations. The MFDs for hydrogen with coupling parameters 
r ee = Fa = 1 and for Al 13+ plasmas with F ee = 0.1 and Fa = 7.19 are shown in Fig. 9 and 10, 
respectively. The dashed curves are the Holtsmark MFDs for a TCP with regularized Coulomb 
potential. Note that the Holtsmark MFD is Z-dependent here (see Eqs. (13) and (24)). To 
demonstrate the importance of attractive interactions we also plotted the MFDs Pq(E) resulting 
from the corresponding electronic and ionic OCPs with F ee and Fa, respectively (open circles). To 
that end the distribution Qo(E) of the total field E = Ei + E2 is calculated as 

Q (E) = j Qe (Ei) Qi (E 2 ) 5 (E — Ei — E 2 ) dEidE 2 = (42) 

from the MFD of the ionic OCP at a charged point Qi (E 2 ) and of the electronic OCP at a neutral 
reference point Q e (Ei). The distribution Qq (E) thus represents the MFD in a TCP assuming that 
the ion-electron attractive interaction is switched off. Here Q e (Ei) and Qi (E 2 ) are taken from 
MD simulations of an OCP. 

Systematic dependencies of the MFD on 5 and F are shown in Figs. 11-16. For fixed F the 
maximum of P(E) shifts only slightly to lower field strengths E with increasing 5, see Figs. 11 and 
12, while the maximum itself increases with 5. This is related to the largest possible single-particle 
field |£^ e (0)| = ei?/2<5 2 , which an electron can produce at the ion. Thus the nearest neighbor 
electronic MFD vanishes for electric fields larger than |£? e (0)|, and smaller 5 will result in larger 
contributions to P(E) at higher fields E with a corresponding reduction of P(E) at small fields. In 
order to demonstrate the enhanced probability of large fields at small 5 and the behavior of PMFEX 
and MD treatments at large fields, the MFD is plotted in Figs. 13 and 14 in a double-logarithmic 
manner. From Fig. 13 it can be deduced that the behavior of the MFD at large fields in H + plasma 



with 5 = 0.2 and T ee = 1 is similar to the nearest neighbor electronic distribution considered in 
detail in Ref. Q. In this case the MFD is strongly reduced at E > \E e (0)\ ~ 12.5E H . 

For fixed 6 and increasing T the MFDs for hydrogen (Fig. 15) and Al 13+ (Fig. 16) show different 
behavior. For hydrogen, like for an ionic OCP, the growing correlations shift the maximum of the 
MFD towards lower electric fields. In the A1 13+ -TCP, P(E) first follows this trend, but then, for 
further increasing T, the maximum turns back to higher field strengths. This can be attributed to 
the growing contribution of the attractive electron-ion interaction and close ion-electron configu- 
rations, which are particularly important for a TCP with highly charged ions. These dependencies 
are very well reproduced by the PMFEX predictions. 

The agreement with the MD data is nearly perfect in most of the studied cases, both for the 
H+-TCP and the A1 13+ -TCP (Figs. 9-16). The PMFEX approximation remains accurate also up 
to high electric fields where the MD data are characterized by strong fluctuations (see Figs. 13 and 
14). Deviations emerge only for strongly coupling situations with large T and a. One example is 
the case of strongly coupled hydrogen with T ee = 1 and 5 = 0.2, i.e. a ~ 6.3 (dotted line and open 
triangles in Fig. 15, see also Fig. 13). Here PMFEX and MD results differ considerably, although 
the HNC treatment is accurate in this case (see Fig. 2). To understand this feature better we recall 
that within PMFEX the Fourier transformed MFD, T(K), fulfills exactly the second moment 
relation (21) in the limit K — > 0. As discussed above only electrons contribute to the second 
moment since the role of ions is negligible (the second term in Eq. (21)). Because of Eq. (36) small 
values of K correspond to large values of the local electric field. One expects therefore that the 
PMFEX yields good results if there are many electrons near the ion. On the other hand, for a 
large electron-electron repulsion T ee = 1 and a light ion like hydrogen, the electrons tend to exert 
only small fields in the ion, for which the quality of the PMFEX is less obvious. 

With increased coupling also the shape of the MFD starts to change. First by a broadening 
of the maximum, and then by the appearance of a shoulder as also reported in Ref. which 
then gets more and more pronounced and finally develops into a second maximum. For the Al 13+ 
TCP in a parameter regime still below the critical values, the MFD is characterized here by the 
formation of the characteristic shape shown in Fig. 17. The HNC approximation still gives the 
correct gu(r) and g e i{f) and the PMFEX well reproduces the broadening and the specific shape of 
the MFD. With respect to g ee (r), however, deviations between HNC and MD emerge (similar to 
those shown for 5 = 0.4a, T ee = 0.2 in Fig. 6), although the electron-electron coupling (T ee = 0.1) 
is still small. The strong ion-electron coupling increases the electron density near the ion which 
introduces additional correlations between electrons, see the discussion above in Sec. VB. This 



will, however, not affect the quality of the PMFEX approximation, since the g ee {r) is not needed 
for the calculation of the MFD at the impurity ion (see Eqs. (13) and (36)). An example for a 
second maximum is given by the strongly coupled hydrogen of Fig. 18, where the parameters are 
close to the critical values. This regime is characterized by the population of bound states and the 
formation of a separate contribution to P{E) at high fields which is mainly due to the electrons. 
Here occur significant deviations between the HNC approach and MD simulations in the RDF 
g e i(r) (see Fig. 3) and the PMFEX approximation cannot predict the shape of the MFD, even 
not qualitatively. But, for coupling parameters, where a classical approach is justified, i.e. when 
bound states are unimportant, the PMFEX approach turns out to be a very reliable method for 
calculating the MFD of a TCP with attractive interaction. 

VI. DISCUSSION AND CONCLUSION 

In this paper our objective was to investigate the microfield distributions in a two-component 
plasmas with attractive electron-ion interactions. Attention has been focused on testing the pre- 
dictions of the PMFEX approximation based on the HNC treatment of static correlations by 
confronting it with the MFDs obtained from MD simulations. One of the basic assumptions of the 
model considered here is the regularization of the attractive Coulomb interaction at short distances 
to introduce quantum diffraction effects in the employed classical approach. 

Two specific rather distinct cases, H + (n e = n«) and Al 13+ (n e = 13nj) two-component plasmas 
with symmetric and largely asymmetric density distributions between plasma species were con- 
sidered. For simplicity we assume bare Coulomb electron-electron and ion-ion interactions while 
the parameter 5 for the regularized ion-electron potential varies from 0.1 to 0.4. The coupling 
strength between plasma particles is measured by the coupling parameters T a p with a, (3 = e,i 
and by the ion-electron potential at the origin in units of ksT, a = T e i/5. Our treatment is 
limited to a parameter regime with a < a c (Z,S), where the critical value a c for 0.1 < 5 < 0.4 
varies 8.32 < a c < 13.4 and 6.66 < a c < 7.33 for H+ and Al 13 + TCPs, respectively. Within this 
parameter regime the g a /3(r) from the HNC equations agree well with the MD simulations. Beyond 
these critical a the HNC equations do either not converge or end up in unphysical solutions while 
the MD simulations remain effective at these strong coupling regimes. A further increase of the 
coupling parameters also leads to the formation of classical strongly bound electronic states with 
no corresponding quantum counterpart. Also the microfield distributions obtained from the HNC 
via the PMFEX approximation agree excellently with the MFDs from the MD simulations except 



of some cases close to the critical T, a. This is somewhat surpri sing since a similar approximation 
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). Therefore we have also 



studied for the OCP deviates from MD simulations (see, e.g., 
tested the PMFEX approximation for an OCP, which in contrast to the TCP turns out to be poor 
when compared with MD simulations although the exact second moment is satisfied within the 
PMFEX. The success of the PMFEX approximation for the TCP is a consequence of the attractive 
interaction and is related to the additional positive electronic part in Eq. (36) which accounts for 
the electric fields created by the electrons at the ions. Obviously, the attractive interactions in 
a TCP favour configurations with large electric fields created at the ion which are well described 
within the PMFEX approximation. On the other hand, in a regime dominated by small local fields 
and hence by small local electronic density the PMFEX deviates from the MD. This feature has 
been clearly observed for a single ion embedded in an electronic OCP in Ref. For the TCP, an 
example is the case of strongly coupled hydrogen with T ee = 1, 5 = 0.2a in Figs. 13 and 15. Here 
some improvement of the PMFEX scheme is required. Such work and the application of PMFEX 
to the case of a neutral radiator are in progress. 
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APPENDIX A: THE EFFECTIVE ELECTRIC FIELDS 



In order to reduce the three-dimensional integration in Eq. (34) to a one-dimensional integra- 
tion and to express the effective fields through scalar potentials u a ji(r) we consider the following 
expression 

J dr l u m {r l )[g a p{\r-r l \)-l]=ATV cf (r, p) [g af} (p) - 1] pdp, (Al) 

where 

Gf(r,p) = ±- [ +P u m {r>ydr>. (A2) 
* r J\r— p\ 

Obviously, the gradient of Eq. (Al) yields the second term in Eq. (34). Consequently, recalling the 
spherical symmetry of the single particle fields, Eq. (34) can be alternatively expressed through 
one-dimensional integrals 

S a (r) = E a (r) + i^L £ qpnp f°° Gf (r, p) [g af3 (p) - 1] pdp. (A3) 
9aR (r) J p 



Here 



Gf(r, P ) = -lGf(r, P ) 



(A4) 



For the regularized Coulomb interaction the last expression yields for r < p and r > p, respectively 



8 PR \fipR/ V'^-/.' 



(A5) 



Gf (r,p) 



(A6) 



For bare Coulomb interaction Eqs. (A5) and (A6) for r < p and r > p are reduced to Go (r, p) = 
and G^P (r,p) = p/r 2 , respectively. 
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FIG. 1: RDFs g a p{r) for a H + plasma with fixed T ee = 0.1 and S = 0.1, 0.2, and 0.4. The lines correspond 
to the HNC approximation while the symbols denote the MD simulations. The different lines and symbols 
represent g|f NC = g™ c (solid lines), g™ c (dashed lines), g^ D = gf^ (filled circles), (open circles). 
The numbers indicate the values of 5. 
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FIG. 2: Same as Fig. 1 with fixed 8 — 0.2 and T ee = 0.01, 0.1, and 1.0. Here the numbers indicate the 
values of T ee . 
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FIG. 4: Same as Fig. 1 for a Al 13+ plasma with fixed T ee = 0.01 and S = 0.1, 0.2, and 0.4. The dotted lines 
and the triangles represent g™ c and g^P, respectively. 




FIG. 5: Same as Fig. 4 for the RDFs g e i(r) and gu(r) with 5 = 0.4, and T ee = 0.01, 0.1, and 0.2 as indicated 
by the numbers. 
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FIG. 6: Same as Fig. 5 for g ee (r). 
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FIG. 7: The second moment (E 2 ) of the MFD (in units E% , sec Eq. (25)) obtained from the HNC scheme 
using Eq. (21) as a function of r ee for a H+-TCP. The dashed line corresponds to the limiting case of an 
ideal plasma (see the text for details). The lines with open and filled circles and squares represent the second 
moments for 5 = 0.1, 5 — 0.2 and S = 0.4, respectively. 
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FIG. 9: Normalized electric microfield distributions for a hydrogen plasma with r ee = Tu = 1 and 5 = 0.4 
as a function of the electric field in units of Eh, Eq. (25). The filled circles represent the MFD from the MD 
simulations and the solid curve the results of the PMFEX. The open circles are the MFD obtained from 
the folding of an electronic and an ionic OCP, see Eq. (42). The Holtsmark distribution (see Eqs. (13) and 
(24)) is shown as a dashed line. 
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FIG. 10: Same as Fig. 9 for a Al 13 + plasma with T ee = 0.1 and T u = 7.2. 
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FIG. 11: Normalized electric microfield distributions for H + -plasmas. The lines with and without symbols 
correspond to MD simulations and PMFEX approximation, respectively. T ee = 0.1 and 5 = 0.1 (solid 
lines), 0.2 (dashed lines), and 0.4 (dotted lines). 






FIG. 14: Same as Fig. 
numbers. 



13 for Al 13 + plasmas with T ee = 0.01 and 5 = 0.1, S = 



0.4 as indicated by the 





FIG. 17: Same as Fig. 12 for 5 = 0.2 and Y ee =0.1. 



